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Introduction 


This report summarizes the research conducted on this contract and reported 
on in several technical reports and papers. Reference 1 is a technical report sum- 
marizing research conducted during Phase I of the project (August 1983 - August 
1985), when research was conducted on the following topics: 

1. Development of a generalized time-domain component mode synthesis tech- 
nique for damped systems; 

2. Development of a frequency- domain component mode synthesis method for 
damped systems, and 

3. Development of a system identification algorithm applicable to general damped 
systems. 

The goals of Phase II of the project (August 1985 - August 1987) were: 

1. To continue development of the multi-input, multi-output system identifica- 
tion method initially studied during Phase I; 

2. To explore the possibility of using multishaker testing to identify parameters 
associated with modeling structures with Lanczos modes; 

3. To continue development of the modal test facility in the Structural Dynamics 
Laboratory by acquiring feedback control hardware; and 

4. To design and fabricate a model suitable for studying the control of flexible 
structures, and to develop associated instrumentation. 

This report presents abstracts of the major publications which have been issued 
on the Phase II topics outlined above. 
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Multi-Input, Multi-Output System Identification 


During Phase I, research was conducted on multi-input, multi-output modal 
testing (or system identification). This work was reported on in Refs. 2 and 3. 
Development of this algorithm was continued during Phase II with the results be- 
ing published in Refs. 4-6. These references provide complete documentation of 
the computer program for State-Space Formulation of Multi-Shaker Modal Analy- 
sis (SFMSMA), including a user’s manual (Ref. 4), fully-documented FORTRAN 
code (Ref. 5), and a paper presented at the 5th International Modal Analysis Con- 
ference (Ref. 6). The principal new features of SFMSMA are: introduction of two 
automatic procedures for reducing model size (order), generalization of the form of 
model damping, and utilization of stable least-squares solution techniques. 

Basic Frequency-Domain Algorithm 

To develop a modal parameter estimation method for general linear, time-invariant 
dynamic systems, it is assumed that the dynamics of the system can be represented 
by a discrete analytical model 

[M]i \t) + [C\m + |A-]*(() = (£>K() (1) 

where [M], [C], and [ K ] are, respectively, the mass, damping, and stiffness matrices 
of order N x jV; [ D ] is an N x N p force distribution matrix; and p(t) represents 
the N p shaker (exciter) forces. x(t) is the vector of N displacements, which is 
assumed to include displacements at all exciter coordinates and all response mea- 
surement coordinates. The order, iV, required to completely describe the dynamics 
of a complex structure is theoretically infinite, but the response of the structure 
at a limited number of measurement coordinates due to excitation over a limited 
frequency bandwidth can be represented by a reduced-order model. The vector x(t ) 
may be “condensed” to a reduced-order vector j(t), and x(t) can be “reconstructed” 
from 7 (t) by the transformations 


7(f) = (W*(i) (2) 

?(f) = hWr(f) + ‘10 - WW (3) 

A reduced-order frequency- domain model is created by applying the Fourier trans- 
form to the time histories in Eq. (1) and then employing Eqs. (2) and (3). The 
final result is 

5M + 1C] (jA 5M + (J?l (^L) fc) = iS)pM (4) 
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Let [af(u;,)] represent a set of N w measured acceleration response vectors at fre- 
quencies u> x, u> 2 , . .., u>n u - Then, reduced acceleration response spectra may be 
obtained by using Eq. (2) to get 

[7(k>i)] = [ 7 M 7 M • • • 7(^jvJ] = [V»c][s(wi)] (5) 

Finally, Eq. (4) may be converted to the following form, from which reduced system 
matrices [C], [AT], and [ D ] can be estimated using a least-squares procedure. 

[ ] 

[[C)[K][D)} [-^(u*)]. =[-fNl (6) 

. [ -PM ] . 

The identified matrices [C] and [K] are then employed in an eigensolution which de- 
termines the eigenvalues (damped natural frequencies and decay rates) and complex 
eigenvectors (complex normal modes). 


Coordinate Reduction and Least-Squares Solution Strategies 

The accuracy achieved in modal parameters for a particular reduced-order model 
hinges upon the selection of [ipc] and [V’n]- Two specific reduction strategies are 
implemented in the SFMSMA computer program: and “independent coordinate 
method” based on Gaussian elimination, and a “principal component method.” 

A widely used reduction technique in finite element analysis and structural dy- 
namics is referred to as Guyan reduction (Ref. 7). It assumes that a subset of 
coordinates may be expressed as a linear combination of the remaining coordinates, 
i.e. if x(t) is partitioned into independent and dependent coordinates 


*> - [ ; 


then the dependent coordinates are given by 

Xd{t) = di]xi{i ) ( 8 ) 

Then, the transformations [0 C ] and [rpp] in Eqs. (2) and (3) become 

hM = |[/)(o)]. M = [ ] (9) 

Equation (8) forms the basis for least-squares estimation of the transformation 
matrix [iftdi] using acceleration spectra. 


[zd(w*)] = hM[x,-(w*)] 
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Equations (9) and (10) thus define the “independent coordinate” procedure for 
model order reduction. 

The second model order reduction strategy provided in the SFMSMA system 
identification program is principal component analysis. Principal component anal- 
ysis seeks transformation matrices [ipc] and [V’h] suc h that the expected value of 
the inner product of the error term in Eq. (3), i.e. 

E{e T (t)e(t)} (11) 

is minimized. In addition to producing the transformation matrices [tpc] and [iPr], 
the approach also provides a measure for use in estimating the appropriate order 
for the reduced-order model. 

Least-squares solutions of linear algebraic equations are required at several 
points in the SFMSMA algorithm, e.g. for determining [C\, [K] and [ D ] from Eq. 
(6) and for determining [ipdi] from Eq. (10). The least-squares solution techniques 
provided are the singular value decomposition (SVD) method and the method of 
Householder transformations. 

Summary 

A modal parameter estimation procedure applicable to linear, time-invariant sys- 
tems has been developed. The technique permits multiple non-coherent input ex- 
citations to be applied to the structure. A consistent set of modal parameters is 
produced by the algorithm, which is formulated such that no assumptions con- 
cerning the form of the linear damping matrix [C] need be made. Even in the 
presence of relatively high noise levels, accurate damped frequencies, decay rates 
and complex modes have been identified for a system with both skew-symmetric 
and non-proportional damping terms. 

Since many practical modal analyses involve numerous measurement locations, 
two different transformation strategies for reducing the size of the model were in- 
vestigated. Both approaches attempt to minimize the amount of required user 
interaction. For the small sample size employed in several example runs, the in- 
dependent coordinate method proved to be more accurate than the principal com- 
ponent method. The latter, however, directly provides an estimate of the order of 
model required to accurately model system dynamics in a given frequency range. 
Two least-squares solution procedures which are noted for their accuracy and sta- 
bility, the singular value decomposition technique and Householder’s method, were 
evaluated. Both methods were found to be accurate and’ stable for this application. 
Due to its faster execution time and relative ease of use, Householder’s method 
appears to be the preferred method. 
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System Identification Based on Unsymmetric 
Block-Lanczos Modeling 

The SFMSMA system identification algorithm described above first identifies 
reduced system damping and stiffness matrices [C\ and [K ] (which are normalized 
so that [M] = I), and these are employed in a standard eigensolution to determine 
system complex modes, frequencies, and damping factors. Several years ago, re- 
search by Wilson (Ref. 8) and Clough (Ref. 9) and their colleagues suggested that 
use of Lanczos vectors, instead of normal modes, could result in more accurate and 
more efficient reduced-order dynamic response calculations. This raised the ques- 
tion of whether Lanczos models would also be useful for other structural dynamics 
applications such as component synthesis and control of flexible structures, and 
whether the parameters associated with the reduced-order Lanczos models could 
be obtained experimentally through some modification of system identification pro- 
grams such as SFMSMA. A companion research project (Contract NAS9-17254 
with NASA Lyndon B. Johnson Space Center) has looked at component synthesis 
and control applications (Refs. 10, 11), while identification of Lanczos models was 
considered under the present MSFC contract. 

Although it has not been possible to complete the development of a system 
identification program based on Lanczos modeling, important preliminary steps 
have been completed. As noted above, the SFMSMA method permits the iden- 
tification of reduced-order modal models of systems with general linear damping. 
To revise this algorithm so that Lanczos models of structures with general linear 
damping may be identified from experimental measurements, it was first necessary 
to extend Lanczos methods for modeling undamped structures (e.g. Refs. 8 and 
9) to accommodate general linear damping. The development of an unsymmet- 
ric block-Lanczos algorithm for application to the dynamics of structures having 
general linear damping, described in Ref. 12, will now be summarized. 

Unsymmetric Block-Lanczos Algorithm 

The second-order equations of motion of the system, Eq. (1), may be converted to 
first order and written in the form 

(A)i(i) + Z(t) = F(t) (12) 

Reference 12 shows that [A] can be reduced to block-tridiagonal form [T], 

r m [g,f i 

[Bi] [M 2 ] [G 2 } t 0 

[F]-'[X][f-] = |T) = • • ■ (13) 

0 [Gr-l] 7 

L [£r-i] m J 
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where 


( 14 ) 

( 15 ) 


1 

i 


[F] = [x u ... t x r ] t [x t }eft nxp 

[F]- t = E = [Y u ...,Y r ], MG» nXp 

[Mi] G ft pXp 

[5,],[Gi] G 3? pXp and upper-triangular matrices 
r : number of blocks 

p : size of blocks 

The columns of [F 1 ] are referred to as right Lanczos vectors , while the columns of 
[E] are referred to as left Lanczos vectors. 

The unsymmetric block-Lanczos algorithm for computing [T] may be summa- 
rized as follows: 

[X,], [n] e given, with [A' 1 ] T (y,| = [/,,] 

[AM = (K,] t [/1][A 1 ] 

For j = 1, . . r — 1 (rp < n) 

[Ml = m - \X,)\Mi\ - [Wj-i][G,-i] T 6 «”*' ([A„][G„] T s 0) 

[PA = - Kmf - K-iI(Bj-if 6 a"" m[a,f = 0) (16) 

[Lj][Uj] = [Pj? [i?j] (L — U decomposition) 

[Bj] s [Uj], [* i+1 ] = [RAM ]- 1 

[Y m ] = 

[Afj+i] = [W[A][X i+1 ] 

Validation Studies 

To test the validity of reduced-order models generated using the unsymmetric block- 
Lanczos algorithm described above, both eigensoltuions and forced-response solu- 
tions were computed for the example 8DOF Beam-Rotor Assembly shown in Fig. 1. 
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Figure 1: 8-DOF Beam- Rotor Assembly. 

Results of these studies are presented in Ref. 12. As a sample of these results, 
Fig. 2a shows the random external force applied at DOF 5, while Fig. 2b shows the 
resulting response at DOF 5. While the original system has 8DOF (16 states), a 
reduced-order Lanczos model with only 4 states accurately describes the response 
if the starting vectors in [JAi] are static deflection vectors based on unit loads at the 
excitation coordinates rather than just arbitrary starting vectors. 



TIME (SEC) 


Figure 2a: External Random Force. 
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TIME (SEC) 


Figure 2.b: Response; Order = 4 (2 x 2). 


Summary 

An unsymmetric block-Lanczos algorithm has been developed and tested. This al- 
gorithm can form the basis for future attempts to use experimental “modal analysis” 
to identify the parameters of reduced-order Lanczos models. The algorithm works 
well for systems having repeated eigenvalues and for systems having general linear 
damping. Computation time can be significantly reduced by use of the Lanczos 
models generated by this algorithm. 
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Control-Structure-Interaction (CSI) Studies 


While the previous two sections of this report have described research consist- 
ing primarily of analytical studies aimed at development of algorithms for system 
identification and reduced-order system modeling, the present section describes the 
development of a control-structure-interaction facility and use of this facility for 
conducting a CSI experiment on a flexible beam test structure. 

Hardware Acquisition 

The principal hardware acquisitions during the course of this project were a Sys- 
tolic Systems PC- 1000 array processor and an Acoustic Power Systems Model 113 
Electro-Seis shaker and Model 114 Dual-Mode power amplifier. Funding for 
these was provided by The University of Texas at Austin. NASA funds originally 
intended for use in upgrading the PC-1000 were transferred to Salaries and Wages 
in order to continue the support of graduate research assistants beyond the original 
contract termination date. This did not adversely impact the research project. 

CSI Hardware and Instrumentation 

Reference 13 describes in detail the development of a control-structure-interaction 
facility including hardware acquisition, design and fabrication of specialized hard- 
ware and instrumentation, and design and fabrication of a flexible beam test struc- 
ture. Figure 3 is a block diagram of the entire CSI system including actuator (A PS- 
113), sensors (accelerometers, velocity transducers, displacement transducers, strain 
gases), system identification hardware GenRad 2515), digital control processor (Sys- 
tolic Systems PC- 1000), and miscellaneous instrumentation and hardware. Figure 
4 shows the configuration of the hanging beam test structure. As indicated on this 
figure, a special cross-flexure was employed to simulate a “pinned-end” condition 
at the top of the beam. However, a detailed finite element model of the flexure was 
required in order to model the support condition accurately. Figure 5 shows the 
APS-113 shaker with displacement and velocity transducers to provide co-located 
feedback signals, an accelerometer for use in cancellation of the mass of the shaker 
armature, and a force sensor mounted on a cross-flexure type stinger. Figure 6 
shows the GenRad 2515, the Systolic Systems PC-1000, and several instrumenta- 
tion “black boxes” used for summing signals and controlling the electromagnetic 
“initial condition” release device. 
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Figure 3: System Block Diagram 
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Figure 4: Test Structure 
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Figure 5: Shaker with Transducers 



Figure 6: GenRad 2515, Systolic Systems PC-1000, and Special Instrumentation 
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CSI Experiment 

The CSI experiment consisted of the following steps: 

1. Fabricate the test structure (hanging beam illustrated in Fig. 4) including 
support flexures and release mechanism; 

2. Attach actuator, sensors, and release mechanism; 

3. Create a detailed finite element model of the beam and cross-flexure support; 

4. Perform an eigensolution to obtain modes and frequencies of the analytical 
model; 

5. Perform an experimental modal analysis to confirm the analytical model; 

6. Develop an optimal feedback controller design based on a single force input 
and multiple outputs including displacement, velocity and strain; 

7. Implement the control using the PC-1000 array processor. 

Table 1 shows the first five natural frequencies obtained from the analytical 
(finite element) model and from the modal test. By including a detailed model of 
the cross-flexure support rather than just assuming a “pinned-end” condition at the 
cross-flexure point, it was possible to obtain agreement to within 2% on the first 
five modes. (Four modes were later used in designing a reduced-order controller.) 


Mode 

Mathematical 

(Hz) 

Experimental 

(Hz) 

Percent 

Difference 

i 

0.311 

0.305 

1.96 

2 

1.796 

1.813 

0.94 

3 

5.402 

5.375 

0.50 

4 

11.030 

10.875 

1.43 

5 

18.605 

18.375 

1.25 


Table 1: Natural Frequencies 


Linear Quadratic Regulator (LQR) design methods were used to determine the 
control gains for the Single-Input, Multiple-Output (SIMO) system. The input u(t) 
is the electrodynamic shaker force, and the measurement (output) vector consists of 
various combinations of the tip displacement and velocity (co-located with input) 
and the strain measurements. The optimal control force is given by 

u(t) = -(B][C]-’f (17) 


13 




where y is the measurement vector 


X T 



€ sl 
e g 2 
e 9 3 
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[C] is determined from the analytical modal responses at the sensors, and [J5] is the 
optimal modal feedback gain matrix derived from LQR theory and computed using 
the PC-MATLAB computer program (Ref. 14). 


Summary 

A control-structure-interaction facility has been developed. The facility includes 
various actuators and sensors, a digital control array processor, and a hanging beam 
test structure. All of the steps in a flexible structure control experiment (analytical 
modeling, verification of analytical model, hardware design and fabrication, calibra- 
tion, etc.) were conducted except for the final step of closed-loop feedback control. 
This step was not completed due to a DC bias problem in the PC-1000 controller. 
Once this problem is corrected the facility should permit many meaningful, CSI 
experiments to be conducted. 
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Overview of Lanczos Modeling In Structural Dynamics 


An overview paper describing several studies of applications of Lanczos methods 
in structural dynamics has been prepared (Ref. 15). Included in this paper is a 
brief description of the unsymmetric block-Lanczos studies described above and in 
Ref. 12. Also described in this paper are methods for component synthesis based 
on Lanczos (or Krylov) modes and for feedback control based on Lanczos modes. 
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Conclusions and Recommendations 


Several significant contributions are contained in the reports and papers pre- 
pared during Phase II of Contract NAS8-35338 and abstracted above. The principal 
contributions are: 

1. The refinement and documentation of a multi-input, multi-output modal 
parameter estimation algorithm which is applicable to general linear, time- 
invariant dynamic systems. 

2. The development and testing of an unsymmetric block-Lanczos algorithm for 
reduced-order modeling of linear systems with arbitrary damping. 

3. The development of a control-structure-interaction (CSI) test facility in the 
Structural Dynamics Laboratory of The University of Texas at Austin. 

Further research on the following topics is suggested: 

1. Combine the unsymmetric block-Lanczos modeling algorithm with the modal 
parameter estimation algorithm to produce a system identification algorithm 
which identifies a reduced-order Lanczos model. 

2. Continue development of the CSI facility and use it to explore methods for 
controlling structures having many closely-spaced-frequency modes and meth- 
ods for controlling multi-component bodies. 
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